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ABSTRACT 

In this paper we posit that galaxy luminosity functions (LFs) come in two fundamentally different types 
depending on whether the luminosity traces galaxy stellar mass or its current star formation rate (SFR). Mass 
function types reflect the older stars and therefore the stellar mass distribution, while SFR function types arise 
from the young stars and hence the distribution of SFRs. Optical and near-infrared LFs are of the mass function 
type, and are well fit by a Schechter function (power law with an exponential cutoff at the bright end). In 
contrast, LFs of the SFR function type are of a different form, one that cannot be adequately described by 
a Schechter function. We demonstrate this difference by generating SFR distributions for mock samples of 
galaxies drawn from a Schechter stellar mass distribution along with established empirical relations between 
the SFR and stellar mass. Compared with the Schechter function, SFR distributions have a shallower decline 
at the bright end, which can be traced to the large intrinsic scatter of SFRs at any given stellar mass. A superior 
description of SFR distributions is given by the "Saunders" function, which combines a power law with a 
Gaussian at the high end. We show that the Schechter-like appearance of UV and Ha LFs, although they are 
LFs of SFR function type, results when luminosities are not corrected for dust, or when average statistical 
corrections are used because individual attenuation measurements are not available. We thus infer that the 
non-Schechter form of the far-IR LFs is a true reflection of the underlying SFR distribution, rather than the 
purported artifact of AGN contamination. 

Subject headings: methods: analytical — methods:numerical — galaxies:evolution — galaxies: fundamental 
parameters — galaxies:luminosity function, mass function 



1. INTRODUCTION 

Schechter (1976) realized that the distribution of opti- 
cal luminosities of cluster galaxies empirically follows the 
same functional form that has been introduced on theoreti- 
cal grounds by Press & Schechter (1974) to describe the halo 
mass function. This functional form is now known as the 
Schechter function. It combines a power law at the faint end 
with an exponential cutoff at the bright end and is uniquely 
determined with three parameters. The function has been 
shown to describe luminosities of galaxies in field environ- 
ments too (Felten 1977). Parameterization of the optical lumi- 
nosity function (LF) using a Schechter function has simplified 
comparisons of different samples of galaxies, including sam- 
ples at different redshifts, and the determination of the cosmic 
luminosity density (e.g., Binggeli et al. 1988). 

With the advent of multiwavelength galaxy surveys, LFs 
began to be constructed in the far-IR (e.g., Lawrence et al. 
1986), near-IR (e.g., Mobasher et al. 1993), UV (e.g., 
Sullivan et al. 2000), as well as for optical emission line lumi- 
nosities (e.g., Gallego et al. 1995). It was generally expected 
that these LFs will also follow the Schechter function, and 
most of these studies, often times using LFs with very lim- 
ited dynamic range, confirmed such expectations. One strik- 
ing exception was the LF in the far-IR (Lawrence et al. 1986; 
Saunders et al. 1990), which was possible to construct over 
a very wide dynamic range (~ 5 dex in space density) and 
which showed a significantly shallower decline at the bright 

1 Visiting Astronomer, Spitzer Science Center, Caltech, Pasadena, CA 
91125. 



end than the exponential decline of the Schechter function. 
While this difference between far-IR and other LFs have been 
acknowledged (e.g., Buat & Burgarella 1998; Takeuchi et al. 
2005), the expectations set by the perceived ubiquitousness of 
the Schechter distribution led some to considered the far-IR 
LF as anomalous and perhaps deviating from the Schechter 
form because of an AGN contamination (e.g., Bothwell et al. 
2011). 

In this paper we show that there is a different explanation 
for such deviations from the Schechter form. We propose 
that there are two fundamentally different galaxy distribution 
functions: (1) of the stellar mass (the mass function, MF 2 ) 
and of the (2) star formation rate (the SFR function). Lu- 
minosity functions in different parts of the spectrum will be 
related more to one type or the other. Optical (especially in 
bands past the 4000A break) and near-IR luminosities arise 
from lower-mass stars that contain most of the stellar mass, 
therefore these luminosity functions belong to mass function 
type. On the other hand the emission in the UV, the nebu- 
lar line emission (e.g., Ha, [Oil]) and the thermal IR is more 
closely related to young stellar populations. Therefore, such 
LFs, if properly dust corrected, should belong to the SFR 
function type. 

What are the true underlying functional forms of the stellar 
mass function and the SFR function? Are they different? In 
the last decade advances in stellar population modeling and 
the availability of large surveys made the determination of the 

2 Since we will only be discussing the stellar mass, we will often be omit- 
ting the adjective "stellar". 
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galaxy stellar masses possible for a large number of galaxies. 
Being more fundamental than the optical luminosity function, 
the MF received prompt attention. It too was found to fol- 
low Schechter's functional form (Cole et al. 2001; Bell et al. 
2003). Indeed it can be said that the LFs in the optical and 
the near-IR reflect the underlying Schechter-like distribution 
of stellar masses. 

Accurate SFR functions are more difficult to construct than 
the MFs due to the caveats and larger uncertainties involved 
in deriving SFRs (e.g., Kennicutt 1998). As pointed out, 
some LFs of SFR type were found to be consistent with a 
Schechter form (UV and Ha) while others, most notably the 
far-IR, was not. The current literature has not fully explained 
this difference. The role of dust has been implied in, e.g., 
Buat & Burgarella 1998; Martin et al. 2005; Takeuchi et al. 
2005; Reddy et al. 2010). A related question, why LFs from 
young stars (UV and Ha) appear to have the same Schechter- 
like distribution as LFs dominated by old stars (optical and 
near-IR), received even less attention. 

This study sets out to determine the intrinsic form of the 
SFR function (in turn testing the adequacy of the Schechter 
function) and to explain why different tracers produce differ- 
ent results. Knowing the appropriate parametric form of the 
SFR function will help interpret the observations at the range 
of redshifts and will facilitate comparison with galaxy forma- 
tion simulations. 

To carry out the search for the functional form of the SFR 
function (SFRF) we adopt a simple framework in which we 
produce volume-complete mock samples of galaxies (§2.1) 
that are described by two quantities: stellar mass and SFR, 
where the stellar mass is drawn from a Schechter function, 
while the SFR is obtained by empirically motivated relations 
between mass and SFR. We apply a series of three such SFR- 
mass relations of increasing complexity, the final of which be- 
ing a relatively realistic representation of the observed SFR- 
mass plane. We then study the SFR distributions produced by 
each relation and discuss functional forms that can be used to 
describe them (§2.2-2.4). Readers not interested in the details 
of these exercises should skip to the summary in §2.5. We 
conclude that Schechter formulation is not adequate for de- 
scribing the SFR function. Instead, functions that replace the 
exponential function at the high end with a Gaussian represent 
a far better description. Next, in §3 we discuss the implica- 
tions of the use of non-Schechter functions for the derivation 
of the SFR density. In §4 we apply the proposed functional 
forms to the observed local SFRF and find excellent agree- 
ment. In §5 we discuss the observed properties of LFs of 
SFR type (UV, Ha and IR) and provide explanation as to why 
UV and Ha LFs appear to be well described using the stan- 
dard Schechter function despite being forms of a SFR func- 
tion. We show that the Schechter-like behavior is a coinci- 
dence stemming from the non-application of dust corrections 
(or from application of only average statistical corrections for 
dust) instead of dust corrections based on more robust mea- 
surements of the attenuation on an individual galaxy basis. 
On the other hand, we show that the non-Schechter form of 
far-IR LFs more closely reflects the true SFR function and is 
not the result of a purported AGN contamination. 

Cosmological parameters Vl m = 0.3, 51a = 0.7, Hq = 
70km s" 1 Mpc -1 are assumed throughout. We express all stel- 
lar masses and SFRs assuming Chabrier IMF. 

2. SFR FUNCTION 



The goal of this paper to evaluate the adequacy of the 
Schechter function for describing the star formation rate func- 
tion and to propose eventual alternatives to this function. In 
order to perform such an evaluation the "true" expected star 
formation rate function must be known. We derive the ex- 
pected SFRF from the combination of two relations which 
are well-determined locally: the stellar mass function, and the 
stellar mass-SFR relation. 

To produce SFR functions we construct mock samples in 
the following way: we draw a large sample from a mass func- 
tion that follows Schechter parameterization. Then, to each 
mock galaxy we associate a SFR based on a SFR-mass rela- 
tion. We explore three different types of stellar SFR-mass re- 
lations: (1) simple power-law relation between SFR and mass 
with no scatter, (2) power-law relation with scatter, and the 
(3) bimodal power-law relation with scatters in both modes. 
As discussed below, each succeeding relation is meant to be 
more realistic than the previous. The last should come very 
close to describing true SFRs. For SFRFs resulting from 
these relations we test the adequacy of the Schechter function 
and search for other functional forms that potentially describe 
them better. 

2.1. Construction of mock samples 

To define the underlying Schechter mass function from 
which the mock samples of galaxies are drawn we adopt pa- 
rameters from Panter et al. (2004) who present a mass func- 
tion based on the SDSS spectroscopic sample. For our 
choice of Hubble constant and IMF these Schechter param- 
eters have the following values: mass function normaliza- 
tion (f>* = 2.7 x 10~ 3 Mpc" 3 dex" 1 , characteristic mass (con- 
verted from Salpeter IMF by dividing the mass by 1.228, 
based on Bruzual & Chariot 2003 models) logM* = 1 1 . 10 (in 
solar mass units, throughout), and the faint-end slope expo- 
nent of a = -1.16. This MF was constructed using data with 
logM* > 7.7, but we verify that its Schechter function fit is in 
excellent agreement with the latest 6 < logM* < 8 MF mea- 
surements from Baldry et al. (201 1), which means that it can 
be safely extrapolated to lower masses. 

From this mass function we draw two volume-complete 
mock samples: (1) a mass-limited sample with M„ > 10 8 M Q 
and (2) a SFR-limited sample with SFR > 0.01 M Q yr _1 . SFRs 
are assigned to each galaxy according to one of the three 
SFR-mass relations as previously described. Distinguishing 
between the two samples with different types of limits is im- 
portant because each affects the shape of the SFR function 
differently. 

The exact choice of limits has no consequence on the in- 
ferences drawn from the analysis, but we wish that they re- 
flect some realistic scenarios. For the mass-limited sample 
we take the limit to be logM* = 8, which is approximately 
the lowest mass for which statistically large samples can be 
extracted from SDSS (e.g., Baldry et al. 2008). To this sam- 
ple we do not impose any limits in terms of SFR. For the 
SFR-limited sample we take the limit of 0.01 M©yr _1 , but 
no mass limit. Throughout this work SFRs represent true, 
dust-corrected SFRs. This SFR limit matches the complete- 
ness limits of surveys that target very nearby galaxies (such 
as the Local Volume Legacy (LVL) survey, Kennicutt et al. 
2008; Dale et al. 2009; Lee et al. 201 1). Most of the galaxies 
that produce stars at the rate around the limit are dwarfs with 
6 < logM* < 8 (Johnson 201 1). This SFR-limited sample will 
probe a fraction of actively star-forming dwarfs at that mass, 
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and principally all star-forming galaxies above logM, = 8 (§ 
2.3). In the local universe this SFR limit probes 99% of the 
total SFR density (§ 3). 

The volume for the mock samples is 10 8 Mpc 3 and was cho- 
sen to be large enough so that the features of the SFR function 
are not significantly affected by the Poisson noise over a large 
dynamic range in space density (~ 5 dex). The mass-limited 
sample contains 3.4 x 10 6 galaxies, and the SFR-limited sam- 
ple up to 7.0 x 10 6 , depending on the SFR-mass relation used. 
Unlike real surveys our samples are volume-complete by con- 
struction, so no completeness corrections are required (for re- 
views of the construction of the observed luminosity functions 
see Johnston 201 1; Takeuchi et al. 2000; Willmer 1997). 

2.2. Relation 1: SFR scales as a power ofM* with no scatter 

Recent studies at low (e.g., Boselli et al. 2001; 
Brinchmann et al. 2004; Salim et al. 2007) and interme- 
diate redshifts (e.g., Noeske et al. 2007; Elbaz et al. 2007) 
have found that for actively star-forming galaxies there 
is a relatively tight and straight sequence in log M* vs. 
log SFR space, which has been dubbed the star-forming 
sequence (Salim et al. 2007) or the galaxy main sequence 
(Noeske et al. 2007). Straight sequence in log space is 
equivalent to a power-law relation between SFR and mass 

(SFR oc M*). The reasons behind the existence of the relation 
are currently the focus of many theoretical studies (e.g., 
Dutton et al. 2010; Dave et al. 201 1). 

Here we adopt the empirical relation derived in the local 
universe (z ~ 0. 1) from Salim et al. (2007) (Eq. 11): 



log SFR = 0.651ogM*-6.33, 



(1) 



where masses and dust-corrected SFRs were obtained through 
the use of UV/optical SED fitting of GALEX and SDSS fluxes 
of galaxies falling in the star forming part of the BPT diagram 
(Baldwin et al. 1981). The relation is shown in Figure 1 (up- 
per panel) and is sub-linear (/? = 0.65). Other studies have 
found different values of /?, but they are usually sub-linear 
(Dutton et al. 2010; Huang et al. 2012). 

Equation 1 is applied to masses drawn from the Schechter 
mass function to obtain their corresponding SFRs. SFRs 
are then binned in 0. 1 dex intervals to obtain SFR functions 
shown in Figure 1 (lower panels). Error bars represent Gaus- 
sian approximation of the Poisson error and are typically ex- 
tremely small for most of the bins, meaning that the features 
of SFRFs are accurately determined. To ensure more even 
weighting when performing the fitting, we add in each bin 
0.03 dex (7%) of systematic error, similar to errors in well- 
determined empirical LFs (Blanton et al. 2001). 

What functional form best describes the SFR functions in 
Figures 1 (lower panels)? Their appearance suggests that they 
would be well fit with a standard Schechter function: 3 



(2) 



where X = SFR, X* is the characteristic SFR, (J)* is the normal- 
ization (expressed in units of Mpc" 3 dex" 1 or Mpc" 3 mag" 1 
throughout) and a is the "faint"-end power-law exponent. We 

3 The logarithmic form of the standard Schechter function is expressed as: 
* s (logA')rf(logX) = ln(lO)0*lO (a+1)<loEX - loEX * ) exp [-10 (lo s x - lo s x *>l rf(logZ). 



show the best fitting Schechter function, obtained by mini- 
mizing x 2 , as solid lines in Figures 1 (lower panels). 4 The 
Schechter fits do not follow the SFRFs exactly. The low-SFR 
("faint" end) slope of both fits tends to be shallower than the 
SFRF points. Similarly, the knee of the fits appears to lie at 
lower SFRs than what is expected visually. Finally, at the high 
end the fits are slightly shallower than the SFR function. This 
mismatch is corroborated with the large \ 2 P er degree of free- 
dom (reduced \ 2 ) values of x 2 = 6.2 and 7.0 for the mass and 
SFR-limited cases, respectively. Why is the Schechter func- 
tion not a perfect fit as may perhaps be expected? 

The answer is that after the power-law transformation, the 
exponential part of the Schechter function becomes modified 
into a Sersic function. Unlike the exponential function which 
has a fixed high-end slope, the Sersic function will have dif- 
ferent slopes based on the extra parameter that is featured in 
it. The reason why the SFR function constructed in this way 
appears to be a Schechter function is because on a logarithmic 
plot the shapes of exponential and Sersic functions are identi- 
cal modulo the scale factor, i.e., we can always pick am scale 
such that the two shapes are exactly the same. 

To properly fit the SFR distribution constructed using Eq. 1 
the Schechter function needs to be modified by introducing an 
additional parameter: the power-law exponent (5 between the 
mass and SFR. We call such function the extended Schechter 
function: 5 



$ ES (X)dX = 



px* 



exp 



<X/X*) 



L//3 



dX, 



(3) 



The exponential part of the standard Schechter function be- 
came the Sersic function, with f3 being equivalent to the Ser- 
sic index. The extended Schechter function is the regular 
Schechter function when (3=1. We confirm that the ex- 
tended Schechter function fits the values of the SFR function 
perfectly and retrieves parameters of the generating MF and 
SFR-mass relation. 6 We do not show these fits in Figure 1 
(lower panels) since they would simply pass through all the 
points with zero deviation. 

The need for extending the Schechter function in order to 
model certain distribution functions has recently been recog- 
nized in Bernardi et al. (2010), following Sheth et al. (2003), 
and in Hopkins et al. (2010). Bernardi et al. (2010) notice that 
galaxy sizes and velocity dispersions scale as power laws with 
respect to the optical luminosity (<fi(L)) and therefore remark 
that "if 4>{L) is well fit by a Schechter function, it makes little 
physical or statistical sense to fit the other observables with 
a Schechter function as well." 7 Hopkins et al. (2010) is the 
only work to our knowledge that has attempted to apply the 
extended Schechter formulation to a SFR-like distribution (a 
simulated IR LF). 

The relation between mass and SFR assumed in this sec- 
tion is very simplistic, so even the extended Schechter func- 

4 Fitting is done in log <j>. 

5 The logarithmic expression for the extended Schechter is: 
*Es(logX)d(logZ) = m(10)|*10 <a ' +1)<1 °e x - lo e x * ) exp [-lu (lo s*- lo e ;<: *>/' 3 ] d(lo 

6 The low-end slope in the extended Schechter function is related to the 
low-end slope of the generating MF, which features in Eqn. 2 as a = B(a' + 
1)-1. 

7 Bernardi et al. (2010) formulation of the extended Schechter function 
(their Eq. 9) appears to have an error (what is listed as \/X should be 1 /X*). 
Furthermore, B in Bernardi et al. (2010) and Hopkins et al. (2010) is the slope 
of mass vs. X and therefore the inverse of our B. 
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Figure 1. Relation 1 and the resulting SFR functions. SFR scales as a sub-linear power of M», with no scatter, the simplest of the three relations that we explore. 
Upper panel shows the dependence of SFR on stellar mass. The relation appears jagged because it is represented as the bivariate density image, which is pixelated. 
The underlying stellar mass function is assumed to have the Schechter form. Dotted lines show limits for mass and SFR-limited samples. Lower panels show 
the resulting SFR functions for the mass-limited (left, logM* > 8), and the SFR-limited samples (right, log SFR > —2). Vertical dashed lines indicate the lower 
cutoff used in the fitting. Red curves represent the best-fitting Schechter functions. Neither Schechter fit describes the SFRF accurately. 



tion will not accurately reproduce all the features of real SFR 
distributions, Nevertheless, it already demonstrates that SFR 
functions cannot be adequately described by Schechter func- 
tions. 

2.3. Relation 2: SFR scales as a power ofM*, with scatter 

While the SFR vs. mass relation is relatively tight, any scat- 
ter around this relation that is not correlated with the mass 
would affect the shape of the resulting SFR function. We 
model this scatter with a Gaussian (in log SFR) of a = 0.4 
dex, again based on the results from Salim et al. (2007). 8 The 
scatter along the SF sequence increases with mass from 0.3 to 

8 Throughout this paper we will refer to Gaussians, keeping in mind that 
in linear SFR these functions are actually log-normal. 



0.4 dex, but the constant value is a reasonable approximation 
for this exercise. This scatter is predominantly intrinsic, since 
the SFR errors in Salim et al. (2007) for galaxies on the SF 
sequence are w 0.2 dex. The SFR-mass relation with scatter 
is shown in Figure 2 (upper panel). 

Figures 2 (lower panels) show the SFR functions for the 
mass and the SFR-limited samples. There are several notable 
differences of these SFRFs with respect to the case with no 
scatter. At the high end the tail is now much less steep ex- 
tending to higher SFRs (note that the horizontal scale in Fig- 
ures 2 (lower panels) is much wider than in Figures 1 (lower 
panels)), so the knee appears "softer". The mass and SFR- 
limited SFRFs start to differ significantly at log SFR < -0.5. 
The SFR function of the mass-limited sample (Figure 2 (lower 
left panel)) now features a turnover and drops off at the low 
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Figure 2. Relation 2 and the resulting SFR functions. SFR scales as the power of M* , with a Gaussian scatter of 0.4 dex in log SFR (log-normal in linear 
SFR). The upper panel shows the dependence of SFR on the stellar mass. Dotted lines show the limits of mass and SFR-limited samples. Lower panels show 
the resulting SFR functions for the mass-limited (left, logM* > 8), and the SFR-limited samples (right, log SFR > —2). Vertical dashed lines indicate the lower 
limits used in fitting. Thick curves (red) represent the best-fitting Schechter functions, which describe the distributions very poorly. Thin curves (blue) are the 
best-fitting extended Schechter functions, which yield very good fits in the fitted regions. 



end. This feature is due to the tail of galaxies close to the 
low-mass limit that scatter below the mean SFR-M* rela- 
tion. The fall off is obviously a mathematical consequence 
of the presence of the mass limit, but it is easy to confuse 
it with volume incompleteness, which is why we show its 
effects separately. For the sample with no mass cut (Fig- 
ure 2 (lower right panel)) there is no such turnover and the 
SFRF continues to rise. Knowing how the sample is selected 
is therefore very important in interpreting the "faint" end of 
any observed SFR function. The fall off at the low SFRs can 
also be seen in the cosmological simulations of the SFR func- 
tion of Dave et al. (2011) and the semi-analytic modeling of 
Fontanot et al. (2012), with both groups using mass limits in 
their simulations. 
In evaluating the functional forms that could be used to de- 



scribe these SFR functions we again start with the regular 
Schechter function (Eq. 2). For the mass-limited sample (Fig- 
ure 2 (lower left panel)) we limit the fitting to the part higher 
than the turnover, (log SFR ^ -0.5, dotted vertical line), since 
obviously the Schechter function (or the extended Schechter 
function) will not be able to reproduce the drop at low values. 
The red line shows the best fit. One can see that the shape 
of the Schechter function is quite inadequate, and the result- 
ing parameters are consequently of little value. The regular 
Schechter function does not perform much better for the case 
of the SFR-limited sample either (red line in Figure 2 (lower 
right panel), with both the low-end slope and the high-end 
drop too steep and the knee too high. Formal reduced \ 2 r val- 
ues are in both cases extremely large (17 and 23, respectively). 
Note however that in some real datasets where the data points 
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have significantly larger error bars and the dynamic range is 
small, a Schechter function fit could be formally acceptable, 
leading one to believe that the distribution is intrinsically of 
Schechter form. 

On the other hand, the extended Schechter fit brings signif- 
icant improvements in describing both samples (blue lines in 
Figures 2 (lower panels)), with = 0.3 and 1.2 for the mass 
and SFR-limited distributions, respectively. Note, however 
that even the extended Schechter function cannot reproduce 
these SFRFs perfectly. The consequence is that the best-fit 
parameters cannot be directly mapped back to parameters of 
the underlying mass distribution and the relationship between 
the mass and the SFR. In SFR-mass relation without scatter, 
/3 represented the slope of the SFR-mass relation. Now, the 
best fits have f3 values of 2.08 and 3.03 for the mass and SFR 
limited distribution respectively, in contrast with SFR-mass 
slope of 0.65. 

The only way to reproduce these SFR functions exactly 
would be to again reverse the process by which the SFRs 
were constructed. This can be achieved with the extended 
Schechter function convolved with a Gaussian. Such a "func- 
tion" would feature 5 or 6 parameters: four of the extended 
Schechter function, the scatter er, and also the mass limit 
for mass-limited samples. We confirm that this function fits 
SFRFs in Figures 2 (lower panels) perfectly (fits not shown), 
with the resulting parameters again having an interpretable 
meaning. However, for any SFRF based on real data the cost 
of two to three additional parameters with respect to the num- 
ber needed to describe the Schechter function will be too large 
- the resulting fits would suffer from a high degree of degen- 
eracy and we will consequently not consider this construct in 
further analysis. 

Bernardi et al. (2010) provided an approximate analytical 
expression for the effect of the measurement error on the ex- 
tended Schechter function in the limit of small a (their Eqn. 
10 and 11). However, that form is not appropriate for the 
level of scatter encountered here, which is dominated by large 
intrinsic scatter. This is exemplified by the fact that the distri- 
bution corrected in that way does not preserve the total num- 
ber density of galaxies since the correction factor is always 
greater than one. 

2.4. Relation 3: SFR is bimodal, each mode scales as a 
power ofM t , with scatter 

After including the scatter in the SFR-mass relationship we 
now add one final element to bring mock SFRs close to the 
realistic ones: galaxy bimodality. Galaxy bimodality is most 
often used to describe the character of optical color distribu- 
tion. The blue mode of the color distribution corresponds to 
the star-forming sequence, which is what was modeled in the 
previous sections. The optical red mode corresponds to galax- 
ies that do not belong to the star-forming sequence and thus 
have little or no SF. We refer to them as passive galaxies. They 
include optically red galaxies with SFRs measurably different 
from zero (e.g., the green valley galaxies detected using UV- 
optical colors; Martin et al. 2007; Salim et al. 2007) and the 
galaxies with upper limits on SFR consistent with no SF. 

If the passive galaxies were taken to have exactly zero SFR 
then the modeling of the SFR distribution reverts to the uni- 
modal case (§2.3), with the only difference being that the 
underlying mass function would be for blue (star-forming) 
galaxies alone. However, it is more realistic to characterize 
passive galaxies with a range of non-zero SFRs, especially 
since passive galaxies on the massive end can reach relatively 



high SFRs (< lMoyr 1 ; Cortese 2012). 

To specify bimodal SFR-mass relations for use in our mod- 
eling we again draw on the data derived in Salim et al. (2007). 
Expressions for separate SF and passive sequences were not 
given in that paper, so we determine them now by fitting 
in each mass bin two Gaussians in log SFR. We find that 
the modeling of the SFR distribution at a given mass with 
two Gaussians describes these distributions remarkably well. 
Unlike in optical color where the colors of passive galaxies 
quickly saturate (the red sequence), the passive sequence is 
quite broad in log SFRs resulting in peaks that are not well 
separated, with no pronounced dip between them. The data 
probe the passive sequence very well by reaching down to 
log(SFR/M„) = -14. The peaks of Gaussians yield the fol- 
lowing SFR-mass relations that are very well described with 
power laws. For the SF sequence: 

logSFR = 0.541ogM*-5.42, (4) 

which is slightly shallower than the SF sequence defined by 
SF galaxies selected using the BPT diagram (Eqn. 1), 
and for the passive sequence: 

logSFR = 0.381ogM„-5.20, (5) 

the scatter of which varies from 1.5 dex at logM* ~ 9 to 0.7 
dex at logM* ~ 11.8, and is again a combination of measure- 
ment errors and, to a larger extent, the intrinsic scatter. For 
simplicity, in our modeling we take the scatter of passive se- 
quence to be fixed at 1 . 1 dex, while for the SF sequence we 
use a scatter of 0.4 dex as previously. 

At each mass we determine the passive fraction from the 
ratio of the area below the Gaussian of the passive sequence 
and the total areas of both Gaussians, which can be obtained 
from: 

~ ^pass^pass 
/pass = ~TZ 77 j \P) 

-WpassCpass + W SF (7 SF 

where N is the height of the peaks. It increases from around 
30% at logM* = 8 to 80% at logM, = 1 1 .5. The passive frac- 
tion is well described as a quadratic function of mass: 

/ pass = 0.0534 log 2 M„- 0.905 logM,+ 4. 144. (7) 

In Figure 3 (upper panel) we show the SFR vs. mass values 
of the simulated bimodal distribution. It was constructed so 
that for each galaxy we first determine if it is passive or ac- 
tive using a random number and Equation 7. Then we draw 
SFRs from appropriate sequence and add the corresponding 
scatter. Since the passive sequence has a fairly large scatter 
we do not allow a passive galaxy to have a SFR greater than 
the SF sequence (Eqn. 1) increased by la (i.e., 0.4 dex). In 
such cases we draw a new value for the SFR. For the mass- 
limited sample a full span in mock SFRs is now ps 9 orders of 
magnitude. 

The resulting SFR functions are shown in Figure 3 (lower 
panels), for the mass and SFR-limited samples respectively. 
The two distributions start to depart below log SFR ps 0. The 
main difference of the bimodal mass-limited SFRF compared 
to its unimodal counterpart is an even more gradual drop to 
lower values due to the relatively low SFRs of the low-mass 
passive galaxies, and the presence of the inflection point (the 
"dip") at log SFR ps — 1, reflecting the bimodal nature of the 
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Figure 3. Relation 3 and the resulting SFR functions. SFRs are bimodal and include the scatter around a star-forming and a passive sequence, with the fraction 
of galaxies in each sequence depending on mass as explained in § 2.4. Upper panel shows the mock SFR-mass relation. SFRs on the passive sequence reach very 
low values. Dotted lines show limits of mass and SFR-limited samples. Lower panels show the resulting SFR functions for the mass-limited (left, logM* > 8), 
and the SFR-limited case (right, log SFR > -2). Vertical dashed lines indicate the lower range used in the fitting the extended Schechter functions (thin, blue 
curves). Extended Schechter functions, as well as Saunders functions (not shown) yield good fits in the fitted regions. For the mass-limited case we also fit a 
double Gaussian (in log SFR, thick green curve), which produces a good fit for the entire SFRF. 



SFR distribution at each mass. On the other hand, the SFR- 
limited distribution remains monotonic because at each SFR 
it is dominated by galaxies in the star-forming sequence. 

We again seek an appropriate functional fit for these SFR 
functions, beginning with the SFR-limited distribution. Since 
the SFR-limited distribution is very similar in the bimodal and 
unimodal case, the regular Schechter function can be ruled out 
as a satisfactory form based on previous considerations. As in 
the case of the unimodal distribution, the extended Schechter 
function represents an excellent fit (Xr = 0.3). We show it in 
Figure 3 (lower right panel) with a blue line. 

While the extended Schechter function represents an excel- 
lent functional form for fitting the SFR-limited SFRF, recall 
that its introduction was motivated by "inverting" the power- 



law mass-SFR dependence in the simple unimodal case with 
no scatter. As a result of that the exponential tail changed into 
a Sersic function. Now, we consider another alternative: a 
function that like Schechter and extended Schechter maintains 
the power law at the low end, but now features a Gaussian (in 
log SFR) at the high end. We refer to it as the Saunders func- 
tion, since it was first proposed in Saunders et al. (1990) to 
model LF at 60 /im. The linear form of the Saunders function 
can be expressed as: 



Note that the "Gaussian" part of this function tends to a Gaus- 
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sian when X > X*, and to a constant when X < X*. This 
modification (the addition of 1 in the argument of log) allows 
the low end to transition smoothly into a power law. 

Fitting the Saunders function to a SFRF constructed using 
the bimodal SFR-mass relation and SFR-limited sample (Fig- 
ure 3 (lower right)) we obtain a very good fit with \ 2 = 0.5. 
This is slightly worse than what the extended Schechter fit 
yielded (\ 2 = 0.3), however, whether one or the other is better 
will depend on the level of scatter, especially in the SF se- 
quence. We know that for the unimodal SFRs, in the limit of 
zero scatter the extended Schechter function is a perfect ana- 
lytic description, but this will be less true for large scatter. For 
the scatter assumed here (0.4 dex for the SF sequence), both 
functions are basically equally good approximations. 

The Saunders function has one significant advantage over 
the extended Schechter function — its parameters are less eas- 
ily perturbed by errors in SFRF and, therefore, the relative 
accuracy of the parameters corresponding to the Saunders fit 
is higher. In the case of the bimodal SFRF discussed here 
the Saunders fit has ~ 2.5 x smaller errors in log SFR* and 
log 4>* and 5x smaller uncertainty in the "faint"-end slope. 
This is the result of the significantly lower level of covari- 
ance among the parameters of the Saunders function. On 
average the Pearson correlation index for Saunders function 
parameters is 0.70, while it is 0.94 for extended Schechter. 
Most importantly, the parameters describing the "faint" and 
the "bright" ends have a correlation of only 0.33 in Saunders 
function (7 and a), and yet 0.87 in the extended Schechter 
function {a' and /?). 

While the Saunders function has been proposed for describ- 
ing the 60 fim LF, which can be considered a type of SFRF, the 
motivation for its introduction was simply to provide a func- 
tional form that better describes the real data than the standard 
Schechter function. Here we show that the reasons behind 
such good representation have to do with the nature of SFR 
and its relation to the stellar mass. 

Turning now to the mass-limited SFRF (Figure 3 (lower left 
panel)), if we only aim to fit the part to the right of the turnover 
(log SFR ^ 0), then the same conclusions hold as in the SFR- 
limited case: the extended Schechter function represents an 
excellent fit (x 2 = 0.5; blue line in Figure 3 (lower left panel)). 
Similarly well does the Saunders function (not shown). Is it 
possible to successfully fit the entire distribution including the 
low-end drop and the dip? Both the extended Schechter func- 
tion and the Saunders function are monotonic and feature a 
power law at the low-end, so they will be incapable to repro- 
duce the inflection or the low-end drop. Thus we test a new 
functional form: a a double (i.e., composite) Gaussian func- 
tion (in log SFR), in which each Gaussian should fit one of 
the two modes of SFR distribution. This function is given in 
the log form as: 



1 



$ GG (logXMlogX) = (<j>* P Gp + (t>* s G s )d(logX), (9) 



where 



1 ( \og 2 {X/X* P ) \ 
G,= ^exp ^-5—1, (10) 



a s V2ir 



exp 



log 2 (X/X*) 
2a\ 



(11) 



represents the Gaussian corresponding to the passive popula- 
tion and 



the Gaussian of the SF population, each with its own standard 
deviation and peak position. Double Gaussian features six 
parameters, but the covariances are weak between each set of 
three. 

Indeed, when we fit the double Gaussian we obtain very 
good results (solid green line in Figure 3 (lower left panel), 
with the individual components shown with dotted green 
lines; on a log-log plot the double Gaussian is represented 
as two parabolas). The fitted Gaussian standard deviations 
are ap = 1.0 and as = 0.5, reflecting the scatters of the pas- 
sive and SF sequences. As in the case of the Saunders func- 
tion, the parameters of the double Gaussian are significantly 
less sensitive to SFRF uncertainties than the parameters of 
the extended Schechter function. While much better than any 
other practical alternative, the double Gaussian is not a perfect 
fit (x? = 1-9). It produces slightly stronger inflection (deeper 
dip) than what is seen in the SFRF. Also, recall that when we 
constructed SFRs we required that SFRs from the passive se- 
quence do not exceed SFRs of the SF sequence by more than 
0.4 dex. However, no such restriction is imposed in the fitting, 
so the passive Gaussian extends a bit too much at the high end. 
Clipping the passive Gaussian where it starts to exceed the SF 
Gaussian brings the reduced x 2 to 1-5- 

Note that the success of the double Gaussian in fitting the 
SFRF is not a mere consequence of the fact that in each mass 
bin we modeled SFR-mass relations as the sum of two Gaus- 
sians because the slope of the SFR-mass relation is signifi- 
cantly larger than zero (i.e., the Gaussians from different mass 
bins have different centers in SFR). 

To our knowledge the double Gaussian was not previously 
considered as a functional form for the SFR function. A sin- 
gle Gaussian (lognormal function in linear SFR) has been 
suggested by Martin et al. (2005). However, such form (a 
parabola on the log-log plot) is apparently inadequate when 
the details of the SFR function are considered, i.e., when it is 
measured precisely below log SFR< 0. Fontanot et al. (2012) 
notice the "double peak" feature in their 0.4 < z < 1.8 mass- 
limited SFR functions and suggest that this "peculiar feature" 
might be connected to bimodality. Our analysis shows that 
bimodality is indeed the explanation. 

2.5. Summary: functional forms describing SFR function 

Summarizing the modeling section we conclude that the 
Schechter function is an inadequate description for any re- 
alistic SFR function. Instead, SFRFs derived from the SFR- 
limited sample can be sufficiently well modeled using either 
the extended Schechter function or the Saunders function (a 
power law with a Gaussian decline at the high end). The lat- 
ter is recommended because its parameters, being less covari- 
ant, will be more stable and yield higher relative accuracy. 
For the mass-limited SFRFs the double Gaussian represents 
an excellent solution. Note that in all of these cases the fitted 
parameters will not directly be interpretable as the parameters 
of the underlying mass function or the SFR-mass relation, 
but they will provide robust descriptions of SFRFs that can be 
compared from one study to another. 

3. STAR FORMATION RATE DENSITY 

There are two primary reasons for which describing the 
SFR function with an analytic function form is useful. One 
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Figure 4. Fraction of the total SFR density that is accounted for by integrat- 
ing mock bimodal SFRF down to a given mass limit (x axis), and SFR limit 
(different curves). To probe 90% of the SFR density in the local universe 
requires sampling galaxies to logM* rj 8 and log SFR > -1. The figure can 
be used to determine the correction factors (1/fraction) to be applied to the 
SFR density determinations in the local universe (z ~ 0.1) obtained from the 
direct numerical integration of the SFRF, and to estimate the yield of planned 
SFR surveys. 

is to characterize this distribution through parameterization, 
to facilitate the study of its evolution with redshift. The sec- 
ond is to use the parameterization to infer, by integration that 
involves extrapolation, the total SFR density. In that case the 
Schechter formulation is especially practical because its cu- 
mulative distribution function is finite when X — > and has a 
simple analytic expression. Since, as we have demonstrated, 
the Schechter formulation does not provide an adequate de- 
scription of the observed SFR function, alternative methods 
are needed for inferring the total SFR density. While all of the 
alternative functions that we have considered (the extended 
Schechter, Saunders and the double Gaussian) also have finite 
cumulative distributions, they cannot be integrated to yield 
functions in a closed form. Therefore, we will instead em- 
ploy our bimodal model (constructed based on z ~ 0.1 SFRs) 
to provide numerical correction coefficients to be applied to 
the SFR density determinations in the local universe obtained 
from the direct numerical integration of the observed SFRF 
down to some SFR limit, based on sample with some mass 
limit. 

The fraction of the total SFR density that will be present 
in the observed local SFRF to a given limit is shown in Fig- 
ure 4. The mass limit can be read continuously from the x- 
axis, while each of the curves represents some SFR limit. The 
black solid curve shows the fraction accounted for at different 
mass limits if there was no SFR limit. Going to logM* = 8 in 
such case (corresponding to the mass limit that we explored 
heretofore) probes 92% of the local SFR density; i.e., the val- 
ues obtained from integrating the SFRF in Figure 3 (lower left 
panel) would have to be corrected upward by 9%. It should be 
noted that in the presence of a mass limit the correction is not 
due to mere extrapolation of the SFRF, but serves to correct 
for lower mass galaxies that contribute at a range of SFRs. 
Having the SFR limit of log SFR = -2 (red dotted curve) leads 
to negligible difference with respect to the case with no SFR 
limit. On the other hand, a limit of log SFR = (green dashed 
curve) captures only 64% of the SFR density in the local uni- 
verse. 



We emphasize that if only the SFR limit is present one could 
recover the missing SFRD by simply extrapolating the faint- 
end power law tail (of course, assuming the slope is well con- 
strained), but if the mass limit is also present one would have 
to account for it too, using, for example, the provided figure or 
performing modeling similar to that presented here. Similar 
technique can also be used to estimate the yield of a planned 
survey or to estimate corrections due to incompleteness. 

4. THE OBSERVED LOCAL SFR FUNCTION 

Next we investigate the shape of the observed SFR func- 
tions and test whether they can indeed be described with the 
functional forms determined based on simulated SFRs in § 
2. We construct the SFR function from UV/optical-based 
SFRs of w 50,000 galaxies from Salimetal. (2007). This 
dataset only has a mass limit (logM* = 8) and can formally 
yield specific SFRs as low as log(SFR/M») = -12, which cor- 
responds to very low SFRs even for massive galaxies. The 
SFRF from constructed from these data is shown in Figure 5 
(upper panel). As expected for the mass-limited sample, the 
SFRF drops off at low values. Overall, its shape is quite sim- 
ilar to the bimodal mock SFRF presented in Figure 3 (lower 
left panel). This is not too surprising because the construc- 
tion of the mock SFRs was guided by these observations, but 
it should be kept in mind that the SFR-mass relations we 
used incorporated some simplifying assumptions (e.g., fixed 
scatter). We previously concluded that the double Gaussian 
function is the most suitable description of the mock bimodal 
SFRF with a mass limit. We fit that function to the observed 
data points and indeed obtain a very good fit (shown by green 
curves in Figure 5 (upper panel)). The resulting goodness of 
fit is x;: = 1 .7. The parameters of the best fit are: 

logSFRp = -2.82, o> = 1.14, log = -2.31 (12) 
log SFR£ = -0.70, cr s = 0.72, log ^ = -1.69. (13) 

To construct the SFRF we co-added 1 /Vmax-weighted prob- 
ability distribution functions for each galaxy's SFR. This is 
possible because the SFRs in Salim et al. (2007) were ob- 
tained using the Bayesian SED fitting. In Bayesian SED fit- 
ting each model SED contributes (proportionally to e~ x I 2 ) to 
the probability distribution of a galaxy parameter (such as the 
SFR). Using full probability distributions for each galaxy's 
SFR (instead of a singular value) has the advantage that it 
produces more realistic distributions for ensemble of galax- 
ies. Since full probability distributions of SFRs are not always 
available, we also considered the SFRF where each galaxy is 
represented by a single value of SFR (the mean of the prob- 
ability distribution). Such SFRF is shown in Figure 5 (lower 
panel). We overlay it with the best double Gaussian fit ob- 
tained using full probability distributions. There is an ex- 
cellent agreement between the two especially in the region 
fit by the star-forming Gaussian, including the high-end tail. 
The latter means that the distribution in the high-end tail is 
not due to some galaxies having broad SFR probability dis- 
tributions reaching very high values. Discrepancies start to 
appear only below log SFR < -3, because galaxies with very 
low SFRs usually have broad (poorly constrained) probability 
distributions, so when these are collapsed into a single value 
for SFR their extent towards very low values gets somewhat 
compressed. In any case, such low SFRs have very little effect 
on any global characterization of the SF. 
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Figure 5. Upper panel. The observed z ~ 0.1 SFR function obtained from 
Salim et al. (2007) data where SFRs where obtained from UV and optical 
SED fitting of an optically selected sample with logM, > 8. SFRF is a com- 
posite of each galaxy's SFR probability distribution, not a single value. Error 
bars are determined from the standard deviation of bootstrap samples. Green 
curves represent double Gaussian fits. Lower panel. Same as Figure above 
except that the SFRF was constructed such that each galaxy's SFR was given 
by a single value (mean of the probability distribution). Green curves are 
repeated from the upper panel. Most of the SFRF stays unchanged. 

Our mock SFR functions in § 2 assumed that masses are 
drawn from a Schechter function. Deviations of the ME from 
the Schechter form are known (Baldry et al. 2008), but they 
are relatively small and mostly pertain to lower-mass passive 
galaxies (Peng et al. 2010). The ability to reproduce the ob- 
served SFRF using functions suggested by mock SFRFs con- 
firms that these deviations are not significant in the context of 
this study. 

5. DISCUSSION 

Our analysis of mock SFRFs and observations from 
Salim et al. (2007) show that real SFRFs have significant de- 
partures with respect to a Schechter function. One usually 
determines the SFR based on the luminosity of a tracer pop- 
ulation of young stars, such as the far-UV continuum lumi- 
nosity, emission line luminosity (e.g., Ha, Oil or Paa), or 
from the dust luminosity in some part of the IR SED (PAH 
lines, mid-IR continuum, far-IR continuum or the total IR lu- 
minosity). It is then to be expected that the luminosity func- 



tions of these various tracers would show similar departures 
from the Schechter distribution as the SFRF. While such de- 
partures have been known for a long time in the IR (especially 
the far-IR; e.g., Lawrence et al. 1986; Saunders et al. 1990; 
Takeuchi et al. 2003), there is a general consensus that the 
far-UV LF as well as the Ha LF do to a large degree agree 
with the Schechter function (e.g., Wyderetal. 2005 for UV 
and Gallego et al. 1995; Ly et al. 201 1 for Ha). How can we 
explain this apparent inconsistency? As the analysis in this 
section will show, the primary reason for this is because the 
observed, uncorrected UV and Ha LFs have a form that is 
similar to a Schechter function by coincidence (similarly for 
LFs where dust is "corrected" by applying average statisti- 
cal relations). On the other hand LFs, where the luminosity 
of each galaxy is individually dust corrected (thus becoming 
a true SFR), do show large departures from Schechter form 
in line with our analysis for SFRFs. The departure of the 
LFs from the Schechter form, even with individually applied 
dust corrections, is more difficult to recognize in high-redshift 
studies which feature small samples with the resulting LFs 
having limited dynamic ranges, which is why this section will 
discuss more robust observational evidence from lower red- 
shifts. 

5.1. UVLFs 

The characterization of galaxies in the ultraviolet has 
greatly improved with the launch of GALEX, which surveyed 
most of the sky in two UV bands (FUV, 1500 A; NUV, 2300 
A). Based on the GALEXI2&F observations of ~ 1000 galax- 
ies Wyder et al. (2005) presented FUV and NUV LFs for local 
(z < 0.1) galaxies. These "early" UV LFs were satisfactorily 
fit with Schechter functions in line with previous UV stud- 
ies. In Figure 6 (upper panel) we present an updated version 
of the FUV LF based on ~30,000 galaxies from Salim et al. 
(2007). All magnitudes are on AB system and are ^-corrected 
to z = and corrected for Galactic reddening. Open points 
show LF of FUV absolute magnitude not corrected for in- 
ternal dust attenuation, i.e., like those of Wyder et al. (2005) 
and many other works that present uncorrected UV LFs. A 
standard Schechter function is fit to the uncorrected LF (red 
curve). While it visually appears as a good fit the large sam- 
ple reveals that \7- is 3.7, a relatively large value. Closer look 
reveals that the faint end-slope of the fit is slightly steeper 
than the LF points and also that the fit falls somewhat more 
steeply at the bright end (M FUV < -20). Wyder et al. (2005) 
were unable to identify these discrepancies because their LF 
had a significantly larger error bars and moreover, because the 
brightest point of their LF was at Mpuv = _ 20, just before the 
departure from the Schechter form starts to become apparent 
at the bright end. Indeed, Schiminovich et al. (2007), using 
the same dataset as the one we use, but analyzed indepen- 
dently, mention the high-end deviation from Schechter. They 
tentatively ascribed it to the AGN contamination in the UV. 
However, broad-line AGNs that could affect the UV contin- 
uum were already removed in these samples, so this expla- 
nation seems unlikely. In any case, this not-so-perfect agree- 
ment between UV LF and Schechter function was generally 
neglected. 

Altogether, the observed UV LF has qualitatively small de- 
partures from the Schechter form, while given our results in § 
2 and 4 we would expect a LF of SFR type (such as the UV 
LF) to be very poorly fit by the Schechter function. This is 
because UV LF, not corrected for dust, is not equivalent to 
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a SFRF. Indeed, the departure from the Schechter form be- 
comes much more severe (x* = 24) when the UV LF is con- 
structed from dust-corrected FUV absolute magnitude, which 
we show as solid dots in Figure 6 (upper panel). The dust at- 
tenuation applied in Figure 6 is obtained on galaxy-by-galaxy 
basis from fitting of the full UV-optical SED (Salim et al. 
2007) and is mostly constrained by the UV slope. If, on 
the other hand one was to apply a fixed dust correction, such 
"corrected" LF would simply be the uncorrected LF shifted to 
brighter magnitudes, while the shape and the steep, Schechter- 
like bright end would remain. Interestingly, even applying 
somewhat more sophisticated statistical dust correction brings 
the "corrected" LF only slightly closer to its true shape. For 
example, many studies apply a mass or optical luminosity de- 
pendent dust correction. From our data on SDSS/GALEX SF 
galaxies we find that the medians in mass bin yield the fol- 
lowing relation between FUV attenuation and stellar mass: 

A FUV = 0.71 logM, -5.16. (14) 

If we apply this statistical correction to FUV absolute magni- 
tudes the resulting LFs still features a relatively steep bright 
end, which in LFs with smaller dynamic range could again 
easily be misinterpreted as conforming to a Schechter func- 
tion. The reason why even the mass-dependent dust correc- 
tion cannot reproduce the true SFRF is because at any given 
mass the dispersion in dust corrections is very large (we find 
0.7 mag scatter in Apuv at logM* = 10.5; similarly large scat- 
ter can be seen Garn & Best 2010 analysis of SDSS Balmer 
decrements; their Figure 4) and it is the values that scatter 
above the average relation that are important in shaping the 
bright end of the LF. 

Guided by our previous considerations regarding the SFRF 
we can attempt to fit more appropriate functions to both the 
uncorrected FUV LF (where the departures from Schechter 
are relatively small, but not negligible) and to the dust- 
corrected FUV LF (where departures are severe). Since the 
GALEX sample is similar to an SFR-limited sample as it is 
Mpuv limited, we try both the extended Schechter function ( 
= power-law + Sersic) and the Saunders function (=power-law 
+ Gaussian). We find that both functions represent good fits 
to both the uncorrected and the dust-corrected FUV LFs, but 
that the Saunders function is somewhat better (x$ of 0.5 and 
1.9 for the uncorrected and corrected LF, respectively, vs. \ 2 r 
of 0.8 and 4.5 for the extended Schechter fit). We show the 
best Saunders fits as green curves in Figure 6 (lower panel). 
Note that the high-end tail of the uncorrected LF is now well 
fit. The parameters of the best Saunders fit are: 

MJ UV = -17.49, cr = 0.31, 7 =-1.23, log<£* = -2.41 

(15) 

for uncorrected FUV LF, and: 

MJ UV = -16.05, cr = 0.73, 7 = -0.83, log(f = -2.24 

(16) 

for the dust-corrected one. 

We now return to the apparent puzzle of why the uncor- 
rected FUV LF is reasonably well fit with a Schechter func- 
tion (left red curve in Figure 6 (upper panel)) when our analy- 
sis has shown that the Schechter function should be an appro- 
priate description for quantities that are proportional to the 
mass and have small scatter with respect to it (for example, 
the optical or the near-IR luminosity), which UV luminosity 
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Figure 6. Upper panel. Far-UV luminosity functions based on Salim et al. 
(2007) sample with (dots) and without (open squares) dust correction. Dust 
corrections are determined on galaxy-by-galaxy basis. Error bars are deter- 
mined from the standard deviation of bootstrap samples. Red lines repre- 
sent best-fitting Schechter functions. Schechter function is an acceptable 
fit only for the uncorrected LF. Lower panel. LFs are repeated from upper 
panel, but the fitting functions (green curves) are now Saunders et al. (1990) 
functions which feature Gaussian high ends instead of exponentials in the 
Schechter function. They provide excellent fits for both the uncorrected and 
the corrected LFs. A LF corrected using average statistical relations (e.g., be- 
tween dust attenuation and stellar mass) would also yield steep Schechter-like 
bright-end slopes, and not the much shallower slope of the true SFRF. 

is not. Here we show that this near match is a coincidence 
arising from two effects that approximately cancel out - sub- 
linearity of the FUV luminosity-mass relation and the scatter 
in that relation. The dust correction is on average larger in 
more massive star-forming galaxies (e.g., Wang & Heckman 
1996; Garn & Best 2010), therefore the slope between FUV 
luminosity not corrected for dust and the stellar mass will 
be even lower (less linear) than between the SFR (i.e., dust- 
corrected FUV luminosity) and mass. Indeed, we determine 
this slope to be j3 = 0.32 (while it was (3 = 0.65 for SFR, i.e. 
dust-corrected FUV luminosity, Eqn. 1). Based on this fact 
alone the Schechter function should be expected to be an even 
worse description for the uncorrected FUV LF. We show this 
in Figure 7. The dashed line shows how the LF would look 
if FUV luminosity was linearly related to mass, without any 
scatter. This is simply a Schechter function. However, the 
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Figure 7. Schematic explanation as to why the Schechter function appears 
to be an adequate description of the uncorrected FUV LF despite the very 
different character of UV and optical/near-IR populations. For a distribution 
to have a Schechter form it needs to have a linear dependence on mass with 
no scatter (dashed curve). Severe sub-linearity, as in FUV luminosity vs. 
mass makes the distribution much steeper (dotted curve), but adding the right 
amount of scatter to such non-linear relation (solid curve) modifies the high- 
end tail into a form that resembles the Schechter function (dashed curve). 
Similar principles would apply to LFs in near UV or the Ha LF. 

actual dependency is sub-linear with j3 = 0.32. Such LF, but 
still with no scatter in L(FUV) vs. mass, is shown as a dotted 
curve. It is much steeper at the bright end than the Schechter 
function, and the high-end tail is shifted to the left. However, 
we also need to take into account the scatter of L(FUV) at a 
given mass. Therefore, we convolve the dotted curve with a 
Gaussian with a = 0.65 mag to obtain the final LF (solid line). 
Convolution has the effect of making the tail shallower again, 
which coincidentally yields a distribution that resembles the 
Schechter function. The differences are relatively subtle: the 
slightly more shallow high-end tail and somewhat less steep 
faint-end slope, and will be easily masked in LFs based on 
up to few thousand galaxies. As we have seen these differ- 
ences can be revealed with LFs constructed from an order of 
magnitude larger samples (Figure 6 (upper panel)). 

This approximate cancelation of sub-linearity and scatter 
happens to hold only for uncorrected luminosity. Since the 
dust-corrected FUV luminosity is less sub-linear with respect 
to the stellar mass the dotted line shifts less to the left. When 
scatter is added to it, it moves the LF more to the right and 
with the shallower slope than that of the Schechter function 
(Figure 8), yielding what can be considered a true SFRF and 
which is well described by Saunders function. 

5.2. Ha and [Oil] emission line LFs 

All of the considerations laid out regarding the UV LFs are 
also applicable to Ha LFs. This is because the Ha lumi- 
nosities also require significant dust corrections to be repre- 
sentative of the true SFRs. Ultimately, the uncorrected Ha 
LF is again only approximately and coincidentally described 
by a Schechter function. Accurate Ha LFs do show de- 
viations. This has been noted by Gilbanketal. (2010) for 
both Ha and [Oil] LFs from SDSS spectra (z ~ 0.1) and by 
Zhu et al. (2009) in their z ~ 1 [Oil] LF from DEEP2 spectra. 
These studies instead fit a double power law (broken line on 
a log-log plot) to the observed emission-line LFs. The dust- 
corrected LFs would show even stronger departures from the 
Schechter form if the dust attenuation is determined for each 



Figure 8. Schematic explanation as to why the Schechter function is not an 
adequate description of the dust-corrected FUV LF (i.e., SFRF). The mean 
relation between luminosity and mass is closer to linear, leading to a smaller 
shift to the left (dotted curve), which results in flatter LF when scatter is 
applied (solid curve) and much larger departure with respect to Schechter 
function (dashed curve). 

galaxy individually. In either case a function such as the Saun- 
ders or the extended Schechter would provide an optimal de- 
scription. 

Paa is emitted in the near-IR so it requires very little dust 
correction (Calzetti et al. 2007). We anticipate that once Paa 
LFs are constructed, they will depart significantly from the 
Schechter form and will feature a Gaussian or a Sersic high- 
end tail instead. 

5.3. IRLFs 

The inadequacy of the Schechter formulation for describ- 
ing IR LFs was noticed already in the mid-80s based on 
IRAS 60 /im data. Thus, Lawrence et al. (1986) proposed 
a double power-law fit (i.e., the exponential cut off of the 
Schechter function was replaced by a less steep power law), 
while Saunders et al. (1990), using an order of magnitude 
larger sample (2800 vs. 300), noticed the curved high-end 
tail of the 60 /Ltm LF so instead recommended the combina- 
tion of a power law for the low end and the Gaussian for 
the high end, which we refer to as the Saunders function 
(Eqn. 8). 9 The Saunders function remains an excellent func- 
tional form for the LFs based on the most recent reductions 
of the IRAS 60 /im data containing five times as many galax- 
ies as used in Saunders et al. (1990), and spanning eight or- 
ders of magnitude in space density (Takeuchi et al. 2003, see 
also Wang & Rowan-Robinson 2010). Despite the success 
of Saunders function some studies continue to use the sim- 
pler double power law (with or without a soft transition be- 
tween the two power laws) to model the far-IR LF, especially 
when smaller sample sizes (< 10 3 ) that probe smaller dy- 
namic range of the LF are involved (e.g., Sanders et al. 2003), 
or when the LFs are constructed for luminosities at A > 60/im, 
where the measurements typically have lower accuracy (e.g., 
Goto et al. 2011). The two forms can be rather similar, and 
distinguishing between them requires very accurate LFs. 

9 A number of papers published since 2005 refer to Saunders function as 
the double exponential, which we find inaccurate and confusing. In math- 
ematics the double exponential refers to functions of the form, f(x) = a b , 
which Saunders function is not. 
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The situation appears to be similar for the mid-IR LFs. 
Shupe et al. (1998) find that the 25 /an LF constructed from 
IRAS data has a shallow high-end tail and can be modeled 
as a double power law. Using ISO 15 /im measurements 
from ISO Xu (2000) and Pozzi et al. (2004) also obtained LFs 
with shallow high-end tails that would be incompatible with a 
Schechter distribution. Given these results it comes as a sur- 
prise that Huang et al. (2007) found that the high end of the 
8/im LF is rather steep and is well fit by a Schechter func- 
tion. They find that this is the case even after they subtract the 
stellar continuum contribution, i.e., when the resulting 8/im 
emission comes primarily from PAH line emission, which is 
considered to be a tracer of young populations and therefore 
of SF (Forster Schreiber et al. 2004; Diaz-Santos et al. 2008). 
Huang et al. (2007) explicitly point out that their results are in 
contrast with Saunders et al. (1990) 60 /tin LF that has excess 
counts above Schechter function. They also point out that 
this agreement with Schechter function cannot be explained 
by their removal of galaxies harboring AGN. Namely, dust- 
obscured AGNs could heat the surrounding dust to relatively 
high temperatures, resulting in the IR SED component that 
peaks in the mid-IR (Fu et al. 2010). Huang et al. (2007) offer 
no explanation for the puzzling inconsistency between their 8 
fjm PAH LF and LFs at longer IR wavelengths. We suggest 
that the reason behind this is because the PAH emission is 
perhaps not a good tracer of the current SF. Namely, there are 
indications that the PAHs exist outside of HII regions where 
they are heated by the general interstellar radiation field pro- 
duced by older stars (Calzetti 2011). By tracing intermedi- 
ate age (~ 1 Gyr) or older stellar populations PAH luminos- 
ity immediately becomes more closely related to the stellar 
mass and less so to the current SF, and thus PAH luminosity 
function can be expected to more closely follow the MF and 
therefore the Schechter distribution. In other words, PAH LF 
is perhaps not a true SFRF 

The curious fact that some LFs of SFR type (primarily UV 
and Ha LFs) are Schechter-like, which we now explain to 
be unrelated to the Schechter form of LFs of mass type, and 
the consequent uncertainties as to the true shape of SFRF 
left a vacuum that some recent studies try to fill by propos- 
ing that the IR LF function is intrinsically also of Schechter 
form, but that some effect not related to SF makes the tail 
assume shallower non-Schechter slope. These studies find a 
culprit among the dust-obscured AGN, which, as mentioned, 
can affect the IR SED, especially in the mid IR. If correct, 
such explanation would disagree with our conclusion (and the 
view of many previous IR studies, e.g., Takeuchi et al. 2003, 
2010; Buat et al. 2007, 2009) that the IR LF, to the extent that 
it measures current SF, and correspondingly the SFRFs, are 
intrinsically not of Schechter form. 

We first discuss the evidence for AGN contamination in 
the mid IR (where it is expected to be stronger) and then 
the far IR (which also dominates the total IR luminosity). 
To test the AGN contamination hypothesis one needs to re- 
move the AGN contribution to the IR luminosities and con- 
struct SF-only IR LF. Fu et al. (2010) used Spitzer IRS spec- 
tra of z ~ 0.7 galaxies to construct LFs at 8 and 15 /im rest- 
frame. Mid-IR spectroscopy allowed them to decompose SF 
and AGN components on a galaxy by galaxy basis (at least 
for the galaxies in the bright tail). AGNs were found to domi- 
nate at the highest 8 and 15 fim luminosities. After removing 
their contribution they find that the resulting LFs are well fit 
by Schechter functions. We note that Fu et al. (2010) probe 
LFs with good precision over only 1 dex of luminosities, so 



the shapes of the resulting LFs are not well determined. Fol- 
lowing in the footsteps of Fu et al. (2010), Wu et al. (201 1) 
use mid-IR spectra to decompose star-forming and AGN con- 
tributions at rest frame 15 and 24 /im for a more local sample 
(z < 0.3). They confirm that after correcting for AGN con- 
tribution the 15 and 24 /im LFs become formally consistent 
with the Schechter functions, however, we notice that their 
data do not probe the high end sufficiently well to rule out 
non-Schechter distribution. 10 Interestingly, Rujopakarn et al. 

(2010) who instead of decomposing the AGN contributions 
completely exclude them find much smaller difference be- 
tween the total and SF-only 24 /im LFs. Consequently, they 
fit both with double power laws and do not consider Schechter 
function at all. To conclude, there is tentative evidence that af- 
ter the AGN correction the mid-IR LF is Schechter-like. This 
possible agreement of the mid-IR LFs with the Schechter dis- 
tribution could be explained if the mid-IR continuum does not 
truly trace the current SFR. Indeed, Kelson & Holden (2010) 
present a model in which most of the mid-IR emission is pro- 
duced by carbon AGB stars with ages between 0.2 and 2 Gyr, 
i.e., while Salim et al. (2009) find that the ~ 15/mi luminos- 
ity is more tightly correlated with the optical B-band lumi- 
nosity, where intermediate age stars dominate, than it is with 
the dust-corrected FUV luminosity of young stars. Significant 
contribution of low mass stars in the heating of mid-IR dust 
would make it be more strongly correlated with the cumula- 
tive SF and therefore the stellar mass, than to the current SF. 
Consequently, the mid-IR LFs may be of mass function type 
and thus be better described by Schechter functions. 

AGN contribution can be accessed relatively directly in the 
mid-IR where the AGN SED peaks, but it is more uncertain in 
the far and the total IR. Wu et al. (2011), extrapolating from 
their 24 /im results, find that the AGN contribution to the to- 
tal IR luminosity is only 10 to 20% (out tologLiR = 11.7), 
yet some studies attempt to correct for AGN contribution in 
the total IR luminosity. This can be done in two ways. One 
is to keep all galaxies but remove fraction of IR luminosity 
believed to come from an AGN, and the other, more crude 
method is to remove galaxies showing signs of AGN from the 
LF altogether. The IR LF constructed using former method 
is equivalent to the SFRF, while the latter will be a lower 
limit to the true SFRF. The recent example of the IR LF that 
completely removes galaxies with AGNs is the one presented 
in Goto et al. (201 1), derived from AKARI data. Goto et al. 

(2011) show full IR LF alongside one that removes AGNs 
based on the number fractions of optically identified AGNs 
in each IR luminosity bin from Yuan et al. (2010). Yuan et al. 
(2010) AGN number fraction is a steeply rising function of 
IR luminosity. Thus the Goto et al. (201 1) LF of non-AGN 
galaxies has a steeper bright end than the full LF, but still not 
so steep to make the authors consider fitting the Schechter 
function instead of the double power law. We confirm this 
by performing the fits ourselves on Goto et al. (2011) LFs. 
Schechter functions yield very poor fit for either the total or 
just the non-AGN LF = 10 and 11 respectively). On the 
other hand, the Saunders functions produce excellent fits with 
X 2 r = 0.6 and 0.7 for the total and non-AGN LFs respectively. 
Since the non-AGN LF is only the lower limit to the real SF IR 
LF, the latter also cannot be a Schechter function. The conclu- 

10 For example, as shown in their Fig. 7 even before the AGN correction 
their 15 fim LF can be fitted by either Schechter or Saunders functions. Un- 
fortunately, one LF point that could help understand what happens at the high 
end was omitted from their SF LF in Fig. 8a. 
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sion that the AGNs are not responsible for the non-Schechter 
form of IR LF is also in line with the exquisite LF at 60 /m 
from IRAS data. Takeuchi et al. (2003) show that the 60 /m 
LF does ha an excess at L(,q > 10 11 5 L , but this excess is only 
ps 0.2 dex in L 6U and lies above the Saunders fit and not just 
above a putative steep Schechter high-end slope. 

Finally, another strong observational confirmation that the 
true SFRF is not of Schechter form comes from the LFs of 
radio galaxies, where the separation between AGNs and the 
star forming component can be achieved more easily than in 
the IR. Thus, for example, Mauch & Sadler (2007) show that 
the 1 .4 GHz LF of star-forming radio galaxies is very well fit 
by a Saunders function, but not Schechter. 

To summarize, while there seems to be some evidence that 
the mid IR LFs after correcting for AGN contribution may 
be consistent with a Schechter function, we propose that this 
could be the consequence of the mid-IR luminosity intrinsi- 
cally tracing SF over longer timescales, and therefore being 
more related to the stellar mass than to the current SF. On 
the other hand, it appears that for either the far-IR or the to- 
tal IR even with very liberal AGN corrections do not produce 
Schechter-like LFs. On the contrary, far and total IR LFs re- 
main well described by functions that replace the exponential 
cutoff of Schechter function with less steep functional forms, 
especially the Gaussian (as featured in Saunders function), in 
line with the expectations from our simulations. 

5.4. Distribution of specific SFRs 

Specific SFR (sSFR) normalizes SFR by stellar mass of a 
galaxy, thus allowing us to asses its SF history and character- 
ize it as bursty, normal or quiescent. The distribution of spe- 
cific SFRs (sSFR function) can offer complementary insights 
to those offered by SFRF alone. We have already shown that 
at any given mass the observed SFR distribution is well de- 
scribed by two Gaussians: one for the SF sequence and an- 
other for a broad, passive sequence. This double Gaussianity 
is preserved when the bimodal SFR-mass relation is projected 
to produce a mass-limited SFR distribution. The same is true 
for the distribution of specific SFRs. It too can be described 
as the composite of two Gaussians in log (SFR/M*). 

Recently, Sargent et al. (2012) have proposed that at z ~ 2 
what we call the SF sequence and model as a single Gaussian 
is in itself composed of two Gaussians - one in which the ma- 
jority of normal SF galaxies lie, and another one, forming a 
bump on the side of the main one, that contains strongly star- 
bursting galaxies. Furthermore, Sargent et al. (2012) present 
a framework in which this two-mode sSFR distribution ex- 
tends to lower redshifts, and they show that it can explain the 
non-Schechter character of the local IR LF. 11 However, our 
analysis of the sSFR distribution of 40,000 z ~ 0.1 galaxies, 
performed in exactly the same way, shows absolutely no indi- 
cation that the SF sequence has a second, starbursting mode. 
Most likely this mode has became negligible since z ~ 2. 
More importantly in the context of this study is that no second 
mode in the SF sequence is needed to produce an IR LF (i.e., 
SFRF) with non-Schechter form. As shown in § 2.3., the non- 
Schechter distribution is primarily the result of a scatter in the 
SFR-mass relation within the unimodal SF sequence alone. 

5.5. Bivariate (s)SFR-M* distributions 

11 To be accurate, Sargent et al. (2012) do not construct a sSFR distribution 
but a sSFR distribution relative to the peak of the SF sequence. The two are 
very similar. 



The approach in this work was to arrive at SFRFs by first 
modeling the SFR vs. M* distribution. SFRF collapses the in- 
formation from the bivariate SFR-M, distribution. Therefore, 
our recommendation is that all studies that report on SFRFs 
or LFs of SFR type should also construct a SFR (or sSFR) 
vs. M* diagram (even a simple scatter plot) to aid in the in- 
terpretation of the SFRF. From such diagram one should try 
to determine the slope and the scatter of the SF sequence and 
get some sense of the fraction of galaxies on the passive se- 
quence, as well as determine the actual mass and SFR lim- 
its. For most purposes even crude mass estimates would suf- 
fice. Techniques developed for the construction of bivariate 
LFs (Takeuchi 2010; Takeuchi et al. 2012; Johnston 201 1) can 
also be applied for the construction of formal bivariate SFR- 
and sSFR-M* distributions. 

6. CONCLUSIONS 

The main conclusions of this study can be summarized as 
the following: 

1. Distributions of the SFR and the stellar mass are fun- 
damentally different. Consequently, the LFs related 
to mass (optical and near-IR) will differ from dust- 
corrected LFs related to SF (UV, emission line, IR). 

2. SFR distributions (SFR functions) are very poorly de- 
scribed by a Schechter functional form, which is ade- 
quate for mass functions. Instead, SFR functions (of 
SFR-limited samples) are very well described by either 
an extended Schechter function (which replaces the ex- 
ponential cutoff with a Sersic function) or by a Saun- 
ders function (where the high end is described by a 
Gaussian in log SFR); see Figures 3 (lower right panel) 
and 6 (lower panel). Both feature four parameters. 
In several empirical cases that we tested (our far-UV 
LF and IR LFs of Goto et al. 201 1 and Takeuchi et al. 
2003) the Saunders function produced somewhat better 
fits. Saunders function has an additional advantage that 
its parameters are less covariant between each other and 
are therefore more robust to LF measurement errors. 

3. SFR functions of mass-limited samples feature a drop 
at the low end even for volume-complete samples. This 
drop should not be confused with incompleteness. The 
shape of SFR functions at the low end is critically sen- 
sitive to the presence of mass limits. Mass-limited SFR 
functions are well described with double Gaussians in 
log SFR (Figures 3 (lower left panel) and 5). 

4. As previous studies have shown, the observed UV 
(and Ha) LFs with no dust correction (or with aver- 
age statistical dust corrections) can approximately be 
described by a Schechter function (Figure 6 (upper 
panel)). Schechter form in them is not fundamental 
(like it is in optical LFs) but is a consequence of two 
effects that by chance approximately cancel each other 
out (Figure 7). Precise LFs (using samples with > 10 4 
galaxies) reveal that even the uncorrected UV LF devi- 
ates from Schechter function. UV and Ha LFs need to 
be dust corrected on galaxy -by-galaxy basis and need to 
have at least a moderate dynamic range for departures 
from Schechter form to be evident. These requirements 
are not always fulfilled in high-redshift studies, lead- 
ing to apparent agreements of UV and Ha LFs with the 
Schechter form. When properly dust corrected, U V LFs 
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follow our mock SFRFs and are successfully described 
by Saunders functions. 

5. LFs in the far and total IR, as well as the radio LF for 
star forming galaxies, behave like our mock SFRFs and 
are poorly fit with a Schechter function even when all 
AGNs are removed. Instead, they are very well fit by 
Saunders functions. 

6. LFs in the mid-lR may intrinsically be Schechter-like 
(after correcting for AGN contribution), which could 
be the consequence of the mid-IR tracing less massive 
(older) stellar populations; i.e., mid-IR LFs are possibly 
not real SFR functions but are more closely related to 
stellar mass functions. Therefore, Schechter functions 
would represent an adequate description. 

7. Whenever possible a bivariate SFR-M* distribution 
should be considered alongside its projections. 

The recognition that optical LFs and MFs can be character- 
ized using a functional form proposed by Schechter (1976) 
has provided extremely valuable guidance in the study of 
galaxy populations across a range of redshifts. As our esti- 
mates of SFRs have improved in recent years, having a simi- 
lar tool to apply to SFR distributions will hopefully add to our 
understanding of galaxies and their evolution. 

S.S.: "I dedicate this paper to the memory of my mother 
Mirjana Makra-Salim (1947-2012). I am forever grateful for 
her selfless love." 
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